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We present the first results of a large spectroscopic sur- 
vey of globular clusters and candidate globular clusters in the 
nearby M 31 galaxy. The survey is aimed at the classifica- 
tion of known candidate M 31 clusters and at the study of 
their kinematic properties. We obtained low-resolution spec- 
troscopy (A/AA a 800 - 1300) for 133 targets, including 76 
yet-to-confirm candidate clusters (i.e. with no previous spec- 
troscopic information), 55 already-confirmed genuine M 31 
clusters, and 2 uncertain candidates. Our observations allowed 
a reliable estimate of the target radial velocity, within a typi- 
cal accuracy of ~ +20 km s . The observed candidates have 
been robustly classified according to their radial velocity and 
shape parameters that allowed us to confidently discriminate 
between point sources and extended objects even from low- 
[ spatial-resolution imagery. In our set of 76 candidates clus- 
£NJ ■ ters we found: 42 newly-confirmed bona-fide M 31 clusters, 
12 background galaxies, 17 foreground Galactic stars, 2 Hn re- 
C$ gions belonging to M 31 and 3 unclassified (possibly M 31 
^) clusters or foreground stars) objects. The classification of a few 

■ other candidates not included in our survey has been also re- 
| assessed on various observational bases. All the sources of ra- 

■ dial velocity estimates for M 3 1 known globular clusters avail- 
| able in the literature have been compared and checked, and a 

■ homogeneous general list has been obtained for 349 confirmed 
OO clusters with radial velocity. Our results suggest that a signifi- 
,— ' cant number of genuine clusters (> 100) is still hidden among 

the plethora of known candidates proposed by various authors. 
Hence our knowledge of the globular cluster system of the M3 1 
galaxy is still far from complete even in terms of simple mem- 
bership. 

^ : 

O ! Key words. Galaxies: individual: M 3 1 - Galaxies:star clusters 
^ - catalogs — Galaxies: Local Group 
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Abstract. 

1. Introduction 

Since the dawn of extragalactic astronomy (Hubble 
the study of the globular cluster (GC) system around the 
Andromeda galaxy (M31) stand out as an important field of in- 
vestigation, providing the ideal counterpart to compare with the 
GC system of our own Galaxy and a testbed for observational 
techniques to be applied to GC systems of more distant galax- 
ies. In present days it is possible to compare directly integrated 
properties and resolved Color Magnitude Diagrams (CMD) 
of M31 clusters, as recently attained with the Hubble Space 
Telescope (HST) (see , e.g. F usi Pecci et al. dl996J) . Brown et 
al. d2004 . Rich et al. d2005l) ). 

In spite of a similar mass and morphology, M3 1 is recog- 
nized to display several differences in its stellar and clusters 
content with respect to the MW. For instance, the halo stel- 
lar populations of the two galaxies (the typical environment 
of most GCs) widely differ in average metallicity; while stars 
in the MW halo are predominantly metal poor (([Fe/H]) 
—1.6 dex (Laird et al. (1988)), those in the outer regions of 
M31 have ([Fe/H]) -0.5 dex ( see, for example, Holland et 
al. dl996h: B ellazzini et al. (2001; Rich et al. i2004 : Durrell 
et al. J2004I) . and references therein). The GC system of M31 
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is much more populous than that of the MW, with more than 
300 confirmed clusters, to compare with the ~ 150 of the MW. 
In particular, Barmby et al. (2001) estimated the total number 
of M31 GCs to be 475 + 25, i.e. more than a factor of three 
larger than in the Milky Way. It is quite clear that understand- 
ing the reasons of such striking differences may shed light on 
the formation histories of the two galaxies and, at large, on the 
general process of galaxy formation. 

Sin ce early systemati c surveys (se e, among others, Vetesnik 
dl962h . van d en Ber gh dl967t Il969l) . Baade & Arp d 1964 . 
Sargent et al. dl977l). Cra mpton et al. dl985t Battistini et al. 
dl980Hl982Hl987Hl993b Sharo v et al. d 19951) . Mochejska et 
al. (|1998), Barmby et al. (|2000j), and references therein), most 
of the known M 3 1 GC candidates (hereafter GCCs) have been 
typically identified by visual inspection of wide-field photo- 
graphic plates, and a more exhaustive analysis with CCD cam- 
eras is still partially lackin g. On t his line, a few recent studie s 
(see, for exampl e, Raci ne dl99lll . Ba rmby & Huc hra d200lh . 
Mo chejsk a et al. d 19981) Hux or et al. d2004l2005k B armbv et 
al. d200lh.P errett et al. d2002h. F usi Pecci et aM2005), Beasley 
et al. d2004 . Puzia et al. (2005)) have clearly shown that this 
kind of new-generation surveys may significantly change our 
knowledge of the number and nature of the clusters harbored 
in the M 3 1 system. To date, th ere are still hundreds of known 
GCCs (-700, see Galleti et al. d2004 . henceforth G04) whose 
nature remains to be ascertained. 

Most of these candidates pertain to the faint-end tail of the 
GC luminosity function whi ch, at present , is far from complete 
(see Fig.[T). Battistini et al. ll98olll982tll987Hl993l) . ranked 
the majority of cluster candidates according to a quality class 
(from "A" to "E" in the sense of decreasing confidence level) 
related to their appearance on the original images. The upper 
panel of Fig.[0shows that the majority of class A and B candi- 
dates have already been confirmed (either spectroscopically or, 
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Fig. 1. Distribution of the original quality-class classifications 
of Battistini et al. (upper panel), and the V luminosity function 
(lower panel) for confirmed M 3 1 clusters (that is all c — 1 
entries in the RBC - see Sect. 1 - according to G04; solid his- 
togram) and candidate M 31 clusters (c = 2 RBC flag; dashed 
histogram). Note the huge number of yet-to-confirm candidates 
at V > 17.0, mainly belonging to class C and D of the Battistini 
et al. classification scale. 



in a few cases, by high resolution images), while most of class 
C and D targets remain to be explored. In any case, one has to 
consider that any comprehensive sample of M 3 1 GCs should 
forcedly rely on coarse databases whose homogeneity is diffi- 
cult to assess. On the other hand, a catalog of M31 clusters as 
complete and homogeneous as possible is clearly overdue, and 
this is a fundamental step for any meaningful comparison with 
other systems. In this framework we have started a long-term 
project to step ahead with this ambitious goal. 

In G04 we re-analyzed the photometric and classification 
data available in the literature, reporting coarse photometry to 
a self-consistent CCD-based magnitude scale, and providing 
infrared information for several hundred GCCs as well, with 
J,H,K magnitudes from the 2MASS database (Skrutskie et al. 
( 2006)). The resulting catalog (Revised Bologna Catalog, here- 
after RBC 1 ) contains at present 1164 entries, including all the 
already confirmed clusters, the known but yet-to-confirm can- 
didates, and the proposed candidates whose "non-M3 1 -GC" 
nature has been definitely established. We believe that keep- 
ing record of the latter objects is very useful all the way, not 
only to avoid duplicating observations, but also to characterize 
the observing properties of the typical contaminants, in order 
to tune up selection criteria and maximize the return of future 
observational campaigns (see G04 for details and discussion). 

For reader's better convenience, it may be useful to recall 
again the RBC classification scheme, as originally used in G04. 



1 The catalog is available 
http : / /www . bo . astro . it/M3 1/ 
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The classification flag c can therefore assume the following val- 
ues: c = 1 confirmed clusters, c = 2 candidate clusters, c = 3 
uncertain candidates, c — 4 background galaxies, c — 5 Hn 
regions, c — 6 foreground stars. In addition, as a result of the 
present discussion (see Sect. 4.3, below), we will also introduce 
here a new class (c = 7) including asterisms/associations. 

As a further line of attack in our work, we are carrying on 
a survey of GCCs located at large (projected) distance from the 
center of M 31, a nearly unexplored realm. A first pilo t run o f 
the survey led to the discovery of B514 (Galleti et al. (|2005), 
hereafter G05), the outermost cluster of the Andromeda galaxy 
ever known. 

Finally we have undertaken a large spectroscopic follow 
up of known candidates, to assess their nature and study their 
physical properties. This is the subject of the present paper, 
where we deal with classification and kinematic of 76 M 31 
GCCs never surveyed before. 

In the following Sec. 2 we briefly introduce the problem of 
M 3 1 GCCs classification, and describe the observational ma- 
terial and its reduction procedures. Sect. 3 is devoted to radial- 
velocity estimates for our sample and the comparison of our 
results with other spectroscopic datasets available in the litera- 
ture. Sect. 4 deals with the classification of the observed candi- 
dates and in Sect. 5 we use the kinematical information for the 
updated sample of confirmed clusters to estimate the mass of 
M 31. Finally, in Sec. 6 we summarize and discuss our results. 

In a companion paper (Galleti et al., in preparation) we 
will complete our study by assessing the problem of a self- 
consistent metallicity scale for M 3 1 GCs, relyin g on ho moge- 
neous measures of the Lick indices (Trager et al. ( 1993)). 

2. Observations and data reduction 

A spectroscopic database of 133 targets has been collected 
through different observing runs along 2004/05, carried out at 
the Roque de los Muchachos (La Palma, Spain) observing fa- 
cilities of the 3.5m Telescopio Nazionale Galileo (TNG) and 
the 4.2m William Herschel telescope (WHT), and at the "G.B. 
Cassini" 1.52m telescope of the Loiano Observatory (Italy). 
Our sample comprises 76 M 3 1 GCCs, while further 55 con- 
firmed clusters (c = 1 RBC flag) were included to consistently 
match our results with other external data sources in the lit- 
erature, together with 2 more questioned objects that need a 
more definitive spectroscopic assessment. The survey mainly 
aims at achieving accurate radial-velocity measurements for a 
wide sample of M 3 1 GCCs through low-resolution long-slit 
and fiber spectroscopy. 

The kinematical piece of information alone can easily dis- 
criminate between genuine M 31 GCs and the most relevant 
class of spurious contaminants, i.e. background galaxies. The 
nature of the latter sources can in fact be recognized because of 
their cosmological recession velocities (typically 6000 < cz < 
50000 km s _1 , see G04), much larger (and opposite in sign) 
than the syste mic vel ocity of M31 (namely V r - -301 km/s, 
Van den Bergh foOOdti ). 

On the other hand, radial velocity alone may not be suf- 
ficient to discriminate between M 3 1 GCs and foreground in- 
terlopers - mostly MW stars - the other major source of con- 
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tamination for our sample. In these cases one has to recur to 
morphological criteria in order to confidently discriminate be- 
tween point (i.e. stars) and extended sources (GCs). To assess 
this important point, we also complemented our spectroscopic 
analysis with supplementary imagery, taken at the Loiano tele- 
scope, for an independent but largely overlapping sample of 86 
RBC objects. We will return on Sec. 4 for a full discussion of 
these results. 

2.1. WYFFOS Data 

Most of our sample has been observed during the nights of 
Nov 21 and 22 2004, using the AutoFib2+Wide Fie ld Fib re 
Optic Spectrograph (WYFFOS; Telting & Corradi fcOOOk 
mounted at the WHT prime focus. We adopted the R1200B 
grating, covering the AA = 4000 — > 5700 A spectral range 
with a FWHM wavelength resolution of 2.2 A. WYFFOS is 
equipped with a mosaic of two EEV-42-80 CCDs, windowed 
and mosaiced such as to have ^ 4300 x 4200 px 2 in total, that 
were read in 2 x 2 binning mode, with a scale of 0.4 A/pixel. 
We used the Small Fibre module, which is made of 150 1 .6 arc- 
sec science fibres, and 10 fiducial bundles for target acquisition 
and guiding. Four different fibre configurations were set up in 
order to target a total of 183 confirmed clusters and GCCs. 

A total of 10 science exposures ranging from 1200 to 
2700 sec were acquired (four in the first night and six in the 
second night). Calibrating exposures included bias, sky and 
lamp flat-field, and He/Ne comparison lamps. For each of the 
four pointings we typically allocated a fraction of ~30% of 
the whole fibre configuration to simultaneously sample the 
sky level. The sky emission was also probed before and af- 
ter each target exposure by dithering the same fibre configu- 
ration. Data were processed in IRAF, using a dedicated pack- 
age (WyffosREDUC) written by Pierre Leisy. The code was 
modified by S.G. to automatically account for the whole data- 
reduction pipe-line; the procedure has been extensively tested, 
comparing the results with those obtained with standard reduc- 
tion techniques for long-slit spectra. The WyffosREDUC pack- 
age is based on the IRAF task DOFIBER, and corrects frames 
for bias, flat-field and for fiber throughputs, extracts spectra for 
each fibre, calibrates them in wavelength, and performs sky 
subtraction as described below. Fibre throughputs have been 
obtained from exposures of sky flat-field frames. 

Sky correction needs special attention when dealing with 
fiber observations; in this regard we explored two independent 
techniques: i) a master sky spectrum was mapped by a grid of 
~35 fibers homogeneously distributed across the whole field of 
view. This output was then subtracted to each science spectrum 
after rescaling for the appropriate fiber throughput; ii) the on- 
target shots were sided (before and after) by two off-target im- 
ages (dithering the telescope 20" E away) to sample sky level 
with the same fiber configuration. No substantial differences 
were found between the two procedures (see Fig.[2]for an illus- 
trative example), and we eventually decided to apply option ( i) 
for our analysis. 

Two M 31 globular clusters (namely B225 and B158), 
for which very accurate radial velocity estimates are available 
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Fig. 2. An illustrative example of the impact of different sky- 
subtraction procedures on WYFFOS spectroscopic observa- 
tions of cluster B039 (upper pane). The raw instrumental spec- 
trum is displayed after sky subtraction i) according to a "master 
spectrum" (Case A) homogeneously sampled across the whole 
field of view by ~ 35 suitably allocated fibres, and ii) by evalu- 
ating sky level from 20" off-target dithered images taken with 
the scientific fiber configuration just before and after the on- 
target shot (Case B). An arbitrary offset in log Flux has been 
added to the data for graphical optimization. 
Lower panel - The flux ratio of Case A and B raw spectra. 
Global standard deviation (per pixel element) amount to a 9% 
(with no evident drift with wavelength), improving to a 4% 
scatter if we restrain to A > 4500 A. 

in the literature (Dubath & Grillmair were observed 

in two different pointings and used as reference templates to 
set radial-velocity zero point. These reference spectra reach 
S /N > 50 per pixel, to be compared with a typical S /N ^ 10 
for the GCC spectra. Two exposures for each scientific target 
have been acquired. The same strategy has been adopted also 
for the runs described below. 

The poor seeing conditions (FWHM ~ 1.5"), in addition to 
thin-cloud sky coverage during both observing nights, and the 
less-than-perfect positioning of some fibres eventually limited 
the final signal-to-noise of several faint targets. For this reason 
only 116 out of the 183 targets ended up with a useful spectrum 
for our analysis. 

2.2. BFOSC Data 

Long-slit spectra for 8 relatively bright (V < 16.5) tar- 
gets were obtained with th e low- resolution spectrograph 
BFOSC (Gualandi & Merighi feOOll) ) operated at the Loiano 
Observatory. The detector was a thinned, back illuminated EEV 
CCD, with 1300 x 1340 px 2 . The observations were carried out 
under typical seeing conditions (FWHM ~ 1.5"), on six nights 
during 2004 (Sep 18-19, Nov 16-17) and 2005 (Jan 02-05). 
With a 1.5 arcsec slit, the adopted setup provided a FWHM 
spectral resolution AA = 4.1 A (A/AA 1300) covering the 
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range AA = 4200 — > 6600 A. We took a He-Ar calibration-lamp 
spectrum after each scientific exposure, maintaining the same 
telescope pointing. Exposure times were typically 45 minutes, 
yielding spectra with a characteristic S fN - 8 per resolution 
element. 

During each observing night and with the same instrumen- 
tal setup, we also collected accurate (S /N > 70) observations 
of four radial-velocity template targets, namely the same M 3 1 
GCs adopted for the WYFFOS run (B158 and B225) plus stars 
HD 12029 and HD 23169 (heliocentric radial velocities from 
the SIMBAD database). Bias, flat field, and sky subtraction 
were carried out using standard packages in IRAF, as described 
in Galleti et al. i 

2.3. DOLORES Data 

The imager/spectrograph DoLoRes at the TNG was used in the 
nights of Sep 8 and Oct 8, 2004, to acquire long-slit spectra of 
9 M 31 GCCs. DoLoRes is equipped with a 2048 x 2048 px 2 
thinned and back-illuminated Loral CCD array providing a 
9.4' x 9.4' field of view. The adopted MRB grism yielded a 
resolution of 6 A FWHM {A/AA = 875) with a 1" slit, across 
a AA = 3800 — » 6800 A spectral range. We typically exposed 
10-15 minutes, reaching S/N - 13 per resolution element. A 
He-lamp spectrum was acquired after each science frame for 
wavelength calibration. Again, during each night we obtained 
good (S/N » 50) spectra for the template clusters B158 and 
B225. The data reduction procedure was the same as for the 
BFOSC data. 



3. Radial velocities 

The heliocentric radial velocities (V r ) of GCCs were ob- 
tained by cross-correlation with the templates spectr a, using 
the IRAF/fxcor package (see Tonry & Davis (1979) for de- 
tails of the technique). We applied a square filter to dampen the 
highest and lowest frequency Fourier components, that heavily 
masked the narrow peaks in the power spectrum. We then fitted 
the power peaks with Gaussians. The typical internal veloc- 
ity errors on a single measure were ~30 km s _1 for WYFFOS, 
~50 km s _1 for BFOSC and ~ 65 km s _1 for DoLoRes spectra. 2 
For each target we cross-correlated two independent spec- 
tra with (at least) four template spectra, obtaining > 8 semi- 
independent estimates of radial velocity. We take the average of 
these values as our final V r , and the standard deviation as our fi- 
nal uncertainty on V r . The typical uncertainties are ~14 km s _1 
for WYFFOS data, -19 km s 1 for BFOSC, and -33 km s 1 
for DoLoRes data. 



2 The standard cross correlation procedure provided poor results for 
background galaxies, given the reduced wavelength range in common 
with the reference template globulars and, in most cases, the presence 
of strong emission lines. For those targets with V r > +6000 km s -1 
we therefore derived the value of cz directly from the measure of a 
few strong emission/absorption spectral features like the [OII]3727, 
[OIII]5007, Ca HK, H/3, H<5, Hy, and MgH lines. This led to sensibly 
higher (but still fully acceptable) uncertainties on the inferred value of 
cz in Table Q 
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Fig. 3. Comparison of radial velocities from the present study 
( Vnhx) with estimates from other authors (H91 = Huch ra et a l. 
( 1991); P02 = Pen-ett etal (l2002t) : F93 = Federici et al. dl993h : 
BOO — Barmby et al. (2000)). In each panel we report the mean 
radial velocity difference and standard deviations between the 
two sets under consideration. 2-cr contours around the mean 
are also reported (dotted lines). 



At the end of the analysis we obtained reliable V r esti- 
mates for 133 targets (116 targets from WYFFOS data, 9 from 
DoLores data, and 8 from BFOSC data), 76 of which were 
previously unconfirmed GCCs. Of the latter, 12 objects were 
eventually classified as background galaxies due to their evi- 
dent cosmological redshift and, in 9 cases, also to striking line 
emission in their spectra (see Sect. 4.2 for further details). In 
addition, two more objects, namely G137 and G270, display 
clear H/3, Hy and [Om] emission lines and a value of V r com- 
patible with M 3 1 Hn regions (see, for a comparison, Diaz et 
al. d!987l) 1. 

The list of the observed targets is reported in Tabled to- 
gether with their apparent V magnitude and V-K color (cols. 2 
and 3, from G04), the estimated V r and the associated uncer- 
tainty (cols. 4 and 5), th e original q uality class (col. 6), accord- 
ing to Battistini et al. ril982l ll987n. and the instrument used to 
obtain the spectra (col. 7). 

3.1. Comparisons with previous studies 

Our subsample of 57 already confirmed M 31 GCs in Tabled 
allows a thorough comparison with other samples of radial ve- 
locities for_M 31 clusters available in the literature [van den 
Bergh ill 9691 henceforth V69); Huchra et al. ( 1982); Huchra 
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Table 1. Radial velocities of observed GCCs and GCs. 
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Fig. 4. Comparisons of radial velocities from low-resolutions 
studies with high-resolution estimates by Peterson (1989) and 
Dubath & Grillmair (1997). The mean differences and standard 
deviations after 2-cr clipping are indicated in each panel. 2-cr 
contours around the mean are also reported (dotted lines). 



et al. dl99ll h enceforth H91); Peters on dl989t) : Dubath & 
Grillm air dl997l) : Federici e t al. dl993l F93); Jablonk a et al 
lll998l J98); Barmby et al. feOOOl B001: Perrett et al. J20Q2i 
P02)]. 

This exercise has the twofold purpose of i) checking the 
reliability of our measures and their overall consistency with 
previous estimates, and ii) probing the mutual self consistency 
of different existing V r datasets and trying to match each of 
them into a common radial-velocity scale. For this we need 
to recover the possible systematic offsets in the V r zero points 
and provide a suitable average for those objects with redun- 
dant/multiple measurements, assessing the intrinsic accuracy 
of the different data sources. 

To do that, it is useful to divide the overall database in two 
categories, i.e. V r measurements coming from low-resolution 
spectroscopy (o-(V r ) - 10-100 km s~'), and those derived 
from high-resolution echelle spectra (cr(V,.) <k 10 km s -1 ). 
In particular, the latter class inclu des th e works of Dubath 
& Grillmair ( 1997) and Peterson ( 1989), who provided ra- 
dial velocities of M 3 1 clusters with a notably high accuracy 
(cr(V r ) < 3 km s~') and excellent agreement for the objects in 
common. We therefore decided to merge both samples into a 
single High Resolution (HR) se t of 24 clusters, referring how- 
ever to the Dubath & Grillmair (1997) V r value for the clusters 
in common. The HR set will provide the backbone of a compre- 



hensive velocity scale joining measures from all the available 
sources. 

In Fig. [5] we present a comparisons of our output for the 
cluster subsamples in common with the low-resolution obser- 
vations by H91, P02, F93 and BOO. In all cases the overall 
agreement is good, within the combined uncertainties of the 
two considered datasets. The comparison with H91 suggests 
the presence of a small systematic, also confirmed by the com- 
parison with HR measures. The agreement with P02 is par- 
ticularly good. Since the typical accuracy of this set is very 
similar to ours, the relatively low dispersion around the mean 
(<x = 20.8 km s~') is a further support to our estimated V r un- 
certainty. The large scatter of the differences with V r estimates 
by BOO is not particularly satisfactory, but the handful of clus- 
ters in common prevents any further conclusion. Unfortunately, 
none of the HR cluster is included in the BOO sample. The com- 
parison of the BOO velocities with the P02 sample (not shown 
here) reveals a good general agreement but also several cases 
of serious inconsistency with 7 out of 21 clusters in common 
exceeding a +100 km s _1 difference in the velocity estimates. 

3.2. Consistency checks: towards a single 
homogeneous list of radial velocities 

The comparisons with the HR set is presented in Fig.0] where 
each panel reports the mean V r difference and standard devia- 
tion for the different samples after a 2-cr clipping iteration. 

The agreement of our V r estimate with the HR ones, for the 
12 clusters in common, is excellent. There is no sizable zero- 
point difference and again the standard deviation witnesses 
the good accuracy of our measures. The same is true for the 
P02 and the F93 sets, although the latter shows a larger scat- 
ter due to a lower intrinsic accuracy of the data. The match 
with the H91 sample confirms, on the contrary, the presence 
of a systematic, as already suggested before and noted by H91 
themselves. We eventually applied the zero-point correction of 
Fig.|4]to this dataset, while we left untouched the P02 and F93 
samples, the zero-point differences being much smaller than the 
involved statistical uncertainties. As for the BOO and J98 data, 
both samples have no clusters in common with the HR sample 
and we had to compare them with other low-resolution obser- 
vations. No systematic offset is needed, although one should 
notice several cases of striking outliers in each sample. 

V69 presented a list of radial velocities for 44 bright M3 1 
GCs and GCCs. The comparison with the HR velocities for the 
19 clusters in common reveals a small systematic (Vv69- Vhr = 
-15.3 +28.0 km s~'), after one outlier rejection. Since all of the 
V69 targets have been re-observed by several authors in more 
recent studies, typically with a better accuracy, we decided to 
retain the V69 estimates only as an external check for contro- 
versial cases. 

Once verified the self-consistency of all the available 
datasets we merged all the sources into one single catalog of 
M 31 cluster radial velocities, tied to the HR set. This re- 
sulted from a weighted average of multiple measures, carefully 
checked on a cluster-by-cluster basis. In general, a weighted 
mean was iterated after 2-cr clipping. When only two incompat- 
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ible measures (at > 2cr) were available, we chose the one with 
lowest uncertainty. The overall compatibility among multiple 
measurements from different sources is, in general, very good. 
We found 18 cases of marginal > 2<x deviation of one measure 
and 9 cases of serious incompatibility, a few of which already 
noted by P02. The latter nine cases probably deserve further 
checks and we shortly mention them: for B104 we adopted 
V r = -395 + 10 km s _1 from BOO, while J98 reports V r = 
+ 120±42kms- 1 ;forB109 we adoptedV, = -372±12kms _1 
from P02 while H91 reports V r = -633.6 + 24 km s" 1 ; for 
B064D 3 we adopted V r = -72 + 10 km s 1 from BOO while 
J98 reports V r = +191.0 + 62 km s _1 ; for B119 we adopted 
the weighted mean between the estimates by P02 and J98, 
V r = -310.1 + 11.1 km s-\ rejecting V r = -137 + 10 km s" 1 
by BOO; for B 124 we adopted V,. = +70 ± 13 km s^ 1 from 
BOO, while J98 reports V r = -75 + 22 km s -1 (the possibility 
of a typo should also be considered in this case); for B301 we 
adopted the weighted mean between the estimates by P02 and 
F93 V,- = -381.9il0.9kms- 1 , rejecting V r = -30+20 kms" 1 
by BOO; for B337 we adopted V r = +50 ± 12 km s 1 from BOO, 
while F93 reports V r = -232 ±26 km s _1 ; for B350 we adopted 
the weighted mean between the estimates by BOO and H91 
V r = -467.3 ± 12.7 km s -1 , rejecting V,- = -251.7 + 26 km s~' 
by F93; finally, for B380 we adopted V r = -13 + 12 km s 1 
from P02 while BOO reports V r = - 1 2 1 ± 3 1 km s" 1 . At least for 
the cases in which only two incompatible estimates are avail- 
able (B104, B109, B124, B337, B380, B064D), a third inde- 
pendent estimate is highly desirable. Also among the 18 cases 
of not-so-strong incompatibility there are four cases in which 
we had to choose between two estimates based on the accuracy 
of the single estimates alone (B008, B047, B144, B314). For 
B131 we obtained V r = -337 ± 3.0 km s" 1 while H91 report 
y r = -444.6 + 28 km s ; for this case we had to recur to V69 
reporting V r = -450 km s _1 , supporting the H91 estimated, 
that we eventually adopted. 

In the following we will use the merged dataset described 
above as our source of V r measures for M 31 GCs. 

4. Classification 

Figure [5] displays the radial velocity distribution of confirmed 
M31 GCs (empty histogram). The distribution is centered on 
the gal axy sy stemic velocity (i.e. V s = -301.0 km s _1 , van den 
Bergh (2000)) and has a standard deviation cr - 160.0 km s . 
The distribution is significantly flatter than a Gaussian curve 
because most of M 31 GCs partake the overall disk rotation 
around the galaxy center (see BOO and references therein). As 
we mentioned in Sec. 2, just on the basis of the kinematical 
piece of information, distant background galaxies can easily be 
excluded as they have typically V r » +300 km s _1 . 

To study the contamination by foreground MW stars we 
obtained a synthetic sample of Galactic stars over the magni- 
tude range 16.0 < V < 19.0, as spanned by our 76 GCCs, in a 
field of 3° x 3° around the po sition o f M3 1 , from the Besangon 
Galactic Model (Robin et al. l2003l ). Down to this magnitude 
limit and across the observed field the synthetic sample pre- 

3 See Updates &Revisions http://www.bo.astro.it/M31/ 
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Fig. 5. Velocity distribution of confirmed M31 clusters. The 
adopted systemic velocity of M3 1 is marked by a dashed seg- 
ment; the dotted segment marks the radial velocity thresh- 
old beyond which the contamination by foreground stars may 
become a serious issue. A Gauss curve with mean [i = 
-301.0 km s _1 and standard deviation cr = 160.0 km s~' is 
also superposed. The dashed histogram is the distribution of 
Galact ic stars as predicted by the Besancon model (Robin et al. 
il2003t) 1. under the assumptions described in Sect. 4. A Gauss 
curve with mean fi = -29.0 km s~' and standard deviation 
cr = 42.6 km s~' is also superposed. 

diets more than 10 4 stars; their radial velocity distribution is 
approximately Gaussian, with mean (V r ) = -29.6 km s~' and 
cr = 42.6 km s , in excellent agreement with the sample of 
confirmed stars in the RBC (G04). 

Since the available sample of confirmed GCs (obviously) 
has not been obtained from the observation of all the sources 
down to V = 19. Oin a 3°x3° field, adirect comparison between 
the two samples would be greatly misleading. As a reason- 
able and conservative normalization, we assume that the global 
catalog of 1164 objects contains as many foreground stars as 
(presently) confirmed GCs having V r estimates, i.e. 313 stars. 
To limit the effects of fluctuations due to low-number statistics, 
we extracted at random from the whole synthetic catalog 100 
samples of 313 stars and we obtained the V r distribution for 
each of them. The shaded histogram is the average of these 100 
distributions. 

From Fig. [5] it can be safely concluded that GCCs with 
V r < -301.0 km s~' cannot be MW stars, hence they must be 
bona fide M31 clusters. According to the above assumptions 
the expected number of MW stars with -301.0 < V r < -150.0 
km/s is < 4, i.e. ~ 1% of the whole sample. On the other 
hand, for V r > -150.0 the contamination by MW stars is 
likely very significant, preventing a fully reliable discrimina- 
tion between foreground stars and M 3 1 GCs based on the ra- 
dial velocity alone. According to these considerations, we de- 
cided to classify as bona fide M3 1 GCs all the candidates with 
V r < -150.0 km s~', requiring further investigation for those 
with higher V r (see Sect. 4.1, below). Note that our approach 
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is more conservative with respect to the majority of previous 
studies that, in general, considered as genuine M 31 GCs all 
the candidates with velocity within ~ +3<x from the systemic 
velocity of M31, in absence of clear evidence of discriminat- 
ing features in the spectra and/or HST imaging revealing the 
stellar nature of the candidate (see, for example, BOO, P02). 
On the contrary, Fig. [5] strongly suggests that GCCs with V r 

around 30.0 + 120 km s _1 (65 objects in the RBC database, 

most of them "confirmed" GCs by virtue of their radial velocity 
alone) should be very carefully considered as they might likely 
include a certain fraction of misclassified MW stars. Clearly, 
any observation assessing the non-pointlike nature of these ob- 
jects would be extremely valuable in this sense (see Sect. 4.1, 
below). 

Using the above described criteria, for our 76 GCCs of 
Table [2 we are left with 

(a) 35 genuine M 31 GCs with V,. < -150.0 km s _1 ; 

(b) 12 background galaxies, with high recession velocity 
and/or line-emission spectra; 

(c) 2 M 31 Hn regions (never observed before); 

(d) 27 candidates, with -150.0 < V r < +100 km s" 1 , possibly 
compatible both with M 3 1 GCs and with MW stars, that 
will be further analyzed in the following. 4 

4.1. Source morphology and foreground star 
contamination 

To further investigate the nature of the 27 GCCs of item "d" 
above we acquired deep white-light BFOSC images of each 
field. In imaging mode, BFOSC has a pixel scale of 0.58"/px 
and a total field of view of 13.0' x 12.6'. All the observations 
have been obtained in 2005, during the nights of Aug 7-9, Sep 
2, 1 1 and 29, Oct 4 and Nov 7. The exposure times ranged be- 
tween 3 and 10 minutes, depending on target brightness and 
atmospheric conditions. The nights were clear (but not photo- 
metric) with a typical seeing around 1.5-2 arcsec FWHM. 

The pointings were accurately chosen in order to include 
in each frame one (or more) of the 27 targets and a number of 
confirmed M 31 GCs. This eventually allowed to extend our 
morphological analysis to a supplementary sample of 56 "con- 
firmed" GCs plus 3 controversial targets. 5 

Images have been bias and flat-field corrected with standard 
reduction procedures. Relative photometry, FWHM and mor- 
phological parameters of each source in the frame - down to a 
5 cr threshold over sky no ise- were obtained with Sextractor 
(Bertin & Arnouts (1996)). Only non-saturated and isolated 
sources were retained in the final catalogs (Sextractor quality 
flag "0"). 

4 According to our constraints on the value of V r , also the contro- 
versial object B341 in TableQ(class c = 3 in the original RBC clas- 
sification) should eventually be comprised in the bona-fide confirmed 
M 31 GCs supporting the original classification by P02. The case of 
B409 deserves, on the contrary, further discussion (see Sec. 14. 1.41 . 

5 We also verified a posteriori that, out of these 59 objects, 13 fur- 
ther entries from TableQwere serendipitously imaged. They are B042, 
B100, B193, B224, B228, B229, B232, B240, B265, B344, B347, 
B366, B045D (see Tabled). 



In spite of the limited spatial resolution of the images, a fair 
assessment of the GCC morphology was made possible by a 
purely differential approach relying on the comparison with ac- 
curate point-spread function (PSF) estimates and on the study 
of the apparent isophotal radius of the detected objects. The rel- 
atively wide field of view of BFOSC allowed in fact a simulta- 
neous high-S/N imagery of the target together with hundreds of 
field stars, and at least a couple of confirmed M 3 1 GCs. Hence, 
the nature of the considered candidate was established by direct 
comparison with surrounding stars and extended sources on the 
same frame. 

In particular, the apparent target size was probed a) in terms 
of its relative excess compared to the local PSF (namely, by 
defining a ratio parameter/? = [target FWHM] /[PSF FWHM]) 
and b) in terms of its departure from the isophotal flux (F) vs. 
isophotal area (A) empirical relationship for stars to be suit- 
ably set in each individual frame according to the overall image 
quality. 6 Both F and A are natural outputs of Sextractor (re- 
spectively FLUXJSO and ISOAREAJMAGE parameters), and 
one could verify empirically that essentially all stars roughly 
obey a A oc F 1 ^ 2 relationship being enclosed within a 0.1 dex 
wide strip in the log F vs. log A plane, corresponding to a vari- 
ation of the measured FWHM of ±5% (see Fig.|6j. Of course, 
hot pixels and cosmic rays appear as off-strip "sub-point-like" 
objects, while any source significantly more extended than the 
stellar PSF stands out in the F-A plot for its larger apparent size 
for a given magnitude level. To some extent, our approach re- 
calls the classical diagnostic tools used, for instance, forhigh- 
redshift galaxy recognition (e.g. Kron ( 1980), Koo et al. ( 1986), 
Molinari, Buzzoni & Chincarini ( 1990)). 

Figure[6]is an illustrative example for a field matching two 
GCCs (open pentagons) and five already confirmed clusters 
(solid triangles). While all the five confirmed GCs are univo- 
cally identified as extended objects in the F-A plane, one of the 
candidates appears as an extended source, while the other one 
is compatible with a stellar point source. Accordingly, from the 
inspection of the individual F-A diagrams, we assigned to all 
the 27 candidates under consideration the flag E (extended) or 
PS (point source). In the following we will refer to this flag as 
the F-A flag or, for brevity, the flag. 

The R parameter and the F-A plane provide two nicely 
complementary tools to judge the candidate extension since 
R is mainly sensitive to the core of the target image, while 
A, as the area within the outermost isophote, is more sensi- 
tive to the wings of the image. This is very well suited to 
study GC candidates since profiles of globular clusters may sig- 
nifican tly var y, depending on their concentration parameter C 
(King ( 1966)). For instance, high concentration clusters may 
have such compact cores that are essentially indistinguishable 
from point sources but a significant difference with respect to 
genuine point sources can be detected looking at the wings of 
the image, sampling the faint halo of the cluster (see, Buonanno 
et al. ( 1981). 



6 Note that the limiting isophote is set at a S /N = 5 ratio per pixel 
element on each frame. This corresponds, in general, to a different sur- 
face brightness magnitude depending on image quality and sky condi- 
tions. 



10 S. Galleti et al.: A large survey of globular clusters in M31. I. Classification and radial velocity. 

Table 2. Morphological analysis for previously confirmed M 31 GCs. 



Name 


Vr ± cr(V r ) 
[kms- 1 ] 


R 


r-A Hag c' ' 


Name 


Vr ± 
[km 


S" 1 ] 


R 


F-A flag 


«(f) 


B213 


-545 


+ 


11 


1.15 


E 1 


B347 


-251 


± 


20 


1.14 


E 


1 


B265 


^t96 


+ 


16 


1.10 


E 1 


B048 


-251 


± 


12 


1.22 


E 


1 


B209 


^60 


+ 


11 


1.18 


E 1 


B220 


-247 


± 


12 


1.75 


E 


1 


B015D 


^45 


+ 


12 


1.31 


E 1 


B167 


-231 


± 


10 


1.15 


E 


1 


B093 


^47 


+ 


12 


1.25 


E 1 


B075 


-212 


± 


12 


1.38 


E 


1 


V031 


^33 


+ 


12 


2.33 


E 1 


B203 


-199 


+ 


12 


1.99 


E 


1 


B161 


—4-13 


+ 


12 


1.20 


E 1 


B154 


-199 


± 


33 


1.40 


E 


1 


B221 


^t06 


+ 


12 


1.24 


E 1 


B188 


-184 


+ 


12 


1.28 


E 


1 


B021 


-403 


± 


12 


1.29 


E 1 


B232 


-182 


+ 


7 


1 .09 


E 


1 


B228 


-457 


± 


26 


1.43 


E 1 


B224 


-161 


± 


2 


2.21 


E 


1 


B100 


-376 


± 


25 


1.45 


E 1 


B200 


-153 


+ 


12 


2.12 


E 


1 


B467 


-342 


+ 


12 


1.26 


E 1 


B184 


-152 


+ 


12 


1.39 


E 


1 


B043D 


-344 


± 


12 


1.51 


E 1 


B367 


-152 


± 


12 


1.07 


E 


1 


B037 


-338 


± 


12 


1.30 


E 1 
















B042 


-338 


± 


10 


1.14 


PS 1 


B103D 


-148 


+ 


12 


1.18 


E 


1 


B401 


-333 


+ 


23 


1.26 


E 1 


B240D 


-148 


+ 


12 


5.83 


E 


1 


B059 


-332 


± 


12 


1.24 


E 1 


B366 


-141 


± 


10 


1.00 


E 


1(3) 


B391 


-325 


+ 


12 


3.73 


E 1 


B272 


-120 


+ 


12 


1.53 


E 


1 


B045D 


-313 


± 


16 


1.22 


E 1 


B355 


-114 


+ 


12 


0.99 


PS 


1(6) 


B024 


-310 


+ 


34 


1 .09 


E 1 


B198 


-105 


± 


12 


1.20 


E 


J 


B222 


-303 


+ 


10 


J ,30 


E 1 


B046 


-98 


+ 


4V 


1 1 3 


r; 




B382 


-302 


± 


12 


1.12 


E 1 


B216 


-93 


± 


10 


1.36 


E 


1(7) 


B164 


-294 


+ 


12 


1.20 


E 1 


B072 


-89 


± 


12 


2.52 


E 




B091 


-290 


± 


12 


1.31 


PS 1 


B190 


-86 


+ 


12 


1.17 


E 




B047 


-291 


± 


12 


1.19 


E 1 


B193 


-59 


± 


2 


1.07 


E 




B354 


-283 


+ 


26 


1.14 


E 1 


B240 


-56 


+ 


5 


1.23 


E 




B210 


-265 


+ 


12 


1.18 


E 1 


B229 


-31 


+ 


5 


1.36 


E 




B214 


-258 


+ 


12 


1.12 


E 1 


B197 


-9 


± 


12 


1.27 


E 




B344 


-252 


± 


13 


1.34 


E 1 

















(t) New classification flag, according to the RBC notation (see also Sect. 1). The newly determined class of the three object for which we modify previous classification is repotted in 
parentheses. 



diagnostics. To this aim we performed a series of tests using 
the half-light radii (r/,) of globular clusters as the characteris- 
tic scale that is more appropriate in this context. According to 
Djorgovski ( 1993), half-light radii of Galactic GCs range from 
0.8 pc to 20 pc, and more than 75% of his sample (118 GC 
with estimates of the half-light radius) have r/, > 2.0 pc. We 
translated these linear radii into angular radii at the distance of 
M31 (D = 783 kpc, see Sect. 5.1, below) and we convolved 
them with the Half Width at Half Maximum (HWHM) of the 
PSF of our images. Finally, we found the minimum /-/, that pro- 
vides a convolved profile larger than that of a point source by 
a factor > 1.05. It turns out that we are able to pick up clus- 
ters having r;, as small as 1.5 pc (corresponding to 0.4", at the 
distance of M31), if the seeing is 2.0" FWHM, and we are sen- 
sitive to r h > 2.0 pc if the seeing is 3.0" FWHM. Note that the 
radii of the isophotes adopted to compute ISOAREAJMAGE 
are typically a factor 2 larger than half-light radii, hence the 
effective sensitivity of the F-A diagram should be significantly 
better than these figures. This is confirmed by a direct com- 
parison that can be made on three M31 clusters included in 
the present programme t hat hav e an estimate of r/, from HST 
data (Barmby & Huchra feOOll) ). namely B167 (r h = 0.33"), 
B232 (r h = 0.66"), and B240 (r h = 0.80"). All of them are 
clearly identified as extended objects in the F-A diagram and 
their R parameters are 1.15, 1.09 and 1.23, respectively (see 
Table 2). Therefore, we conclude that the image analysis tech- 
nique adopted here is able to recognize the extended nature of 
the large majority of M3 1 clusters even with low -resolution im- 
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Fig. 6. F-A diagram for one of the images we analyzed (as 
an example). All the detected sources are plotted as small 
points. Filled triangles correspond to already confirmed clus- 
ters, empty pentagons to the candidate clusters under consider- 
ation. The continuous lines enclose the locus of points sources. 
Extended objects are expected to populate the upper left corner 
above the strip, while cosmic spikes and CCD cool/hot pixels 
must be confined to the lower right corner of the plot. 

It may be useful to have an idea of the sensitivity of the 
adopted technique, i.e. what is the size of the most compact 
M31 GC that can be recognized as an extended object by our 
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ages, provided that the target image is of high signal-to-noise 
so that the image profile is well constrained. 

4.1.1. Checking previously confirmed "genuine" 
clusters with V, < -150.0 

Table |2 reports the results of the image analysis described 
above for the M 3 1 clusters previously confirmed by other au- 
thors that were included in our BFOSC imaging survey. The 
first 42 entries in the table concern clusters that have V r < 
-150.0 km s _1 , and can be unambiguously identified as gen- 
uine GCs on the basis of their radial velocity alone, accord- 
ing to our previous arguments. We note that all of them have 
R > 1.07, and all except three have R > 1.10. Only two clusters 
obtain the PS flag based on their position in the F-A diagram, 
but both of them have quite large R, 1.14 and 1.31 for B042 
and B091, respectively. Hence our morphologic criteria cor- 
rectly recognize all the genuine M 3 1 GCs considered here as 
extended objects. 

Also based on these results, we can complete our classifi- 
cation scheme and devise the following supplementary criteria 
for GCCs with V r > -150.0 km s 1 : 

(a) a candidate is classified as a "bona fide M31 cluster" (c = 1 
in Tables 2, 3, and 4) if R > 1.10, independently from the 
assigned flag, or if R > 1 .05 and F-A flag "E". 

(b) it is classified as a "bona fide star" (c = 6) if R < 1.05 and 
F-A flag "PS" orifR< 1.0. 

(c) it is considered as an "uncertain object" (c = 3) if 1.0 < R < 
1.05 and flag "E", or 1.05 < R < 1.10 and flag "PS". These 
non-decidable cases can be resolved only with further data 
and/or analysis. 

The overall classification picture for the whole GCC sam- 
ple of Table 2 and the following ones is summarized in Fig.0 

4.1 .2. Image analysis for previously "confirmed" 
clusters with V r > -150.0 km s _1 

The last 14 entries of Tab. [2] concern clusters previously con- 
firmed by other authors, that have V,- > -150.0 km s _1 , i.e. a 
range of velocity that may suffer from contamination by MW 
stars. As noted in Sect. 4., those that have been confirmed only 
by virtue of their radial velocity may deserve further analysis 
to obtain a firm classification. According to the above criteria, 
12 of these candidates are classified as genuine M31 GC, one 
(B366) is classified as "uncertain" and one (B355) is classified 
as "bona fide star". It is interesting to note that B355 was clas- 
sified as a confirmed GC by Perrett et. al. (2002) relying only 
on its radial velocity. Finally, B216 makes a case on its own 
that will be discussed in detail in Sect. 4.3. 

4.1 .3. Image analysis for candidates with 
V r > -150.0 km s" 1 

The results of the image analysis for the 27 GCCs from our 
sample of Tabled with V r > -150 km s are reported in 
Table |3] According to the above criteria, 7 of them are classi- 



Table 3. Morphological analysis for the 27 M31 globular clus- 
ters candidates of Table^with V r > -150.0 km s . 



Name 


Vr ± cr(V r ) 
[kms- 1 ] 


class' 01 


R u» 


F-A flag ,rl 


C WJ 


B189 


-148 


± 


8 


A 


1. 26 


E 


1 


B162 


-146 


+ 


8 


A 


1.27 


E 


1 


B217D 


-136 


± 


9 


D 


0.94 


PS 


6 


B187 


-130 


± 


10 


A 


1.47 


E 


1 


B371 


-127 


± 


18 


B 


1.38 


E 


1 


B255D 


-107 


± 


14 


D 


2.5 


E 


1 


B126D 


-92 


± 


18 


D 


1.00 


PS 


6 


B090D 


-94 


± 


8 


D 


1.15 


E 


1 


B260D 


-93 


± 


6 


D 


1.03 


PS 


6 


B362 


-81 


± 


4 


A 


1.09 


PS 


3 


B330D 


-62 


± 


18 


D 


1.01 


PS 


6 


B036D 


-54 


+ 


7 


D 


0.99 


PS 


6 


B243D 


-53 


± 


18 


D 


1.03 


PS 


6 


B388 


-50 


± 


9 


B 


1.34 


E 


1 


B027D 


-50 


+ 


7 


D 


0.99 


PS 


6 


B221D 


-50 


± 


8 


D 


1.05 


PS 


6 


H126 


-19 


± 


10 


/ 


1.06 


PS 


3 


NB70 


-17 


± 


11 


E 


1.02 


PS 


6 


B275D 


-13 


+ 


6 


D 


0.98 


PS 


ft 


B073D 


-12 


± 


13 


D 


1.02 


PS 


6 


B226D 


-9 


+ 


8 


D 


1.00 


PS 


6 


B109D 


-10 


+ 


15 


D 


1.03 


PS 


6 


B021D 


6 


± 


7 


D 


0.95 


PS 


6 


NB65 


g 


+ 


6 


E 


1.03 


PS 


6 


B237D 


10 


+ 


8 


D 


1.03 


PS 


6 


B473 


1 1 


+ 


7 


E 


1.00 


PS 


6 


B213D 


19 


± 


9 


D 


1.06 


PS 


3 



(a) Quality flag according to the Battistini et al. |l980; 1982; 198^ [T99l classification 
scale; 

(b) Target apparent size, relative to the PSF (namely R = [target FWHM]/[PSF FWHM]; 

(c) E(xtended) or P(oint) S(ource) classification obtained from the F-A diagram; 

(d) RBC classification flag as in Tab.[5](see Sect. 1.). 



Table 4. Candidates with controversial classifications 



Name 


Vr ± 
[km 


o-(K) 


R 


F-A flag 


c 


B055 


-308 + 


8 


1.27 


E 


1 


B121 


-24 + 


33 


1.08 


PS 


3 


B409 


-40 + 


10 


1.00 


PS 


6 



fied as genuine GCs (namely B162, B187, B189, B371, B388, 
B090D, and B255D), 17 are flagged as stars and 3 remain un- 
certain. We note that all the newly confirmed clusters are very 
clearly recognized as extended: all of them have flag=E and R 
ranges between 1.15 and 2.5. Note also that 5 of the 6 candi- 
dates in the range -150 < V r < -100 km s _! turned out to be 
genuine GCs, while for the 21 targets with V r > -100 km s _1 
we only detect 2 clusters and 16 bona fide stars. This fur- 
ther confirms how severely star contamination can affect in the 
range of V r around 30 + 130 km s _1 . 

4.1 .4. A few controversial objects 

The three entries of Table |4] concern candidates that obtained 
different classifications from different authors. B055 was clas- 
sified GC by P02, based on its radial velocity, while BOO clas- 
sified it as a star. Our image analysis indicates unambiguously 
that B055 is an extended object, thus confirming the classifica- 
tion by P02. B121 was classified GC by H91 and "star" by BOO 
and our analysis does not help to resolve the controversy. 

Object B409 was classified as a GC by F93, based on its 
radial velocity alone, and it was classified as a background 
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Fig. 7. Summary of the adopted classification criteria illustrated in the R vs. V r diagram, where different symbols are used 
according to the value of the F-A flag of the targets. Circles represent extended objects (F-A flag = E) while triangles are point- 
like sources (F-A flag = PS). Bona fide stars (c=6) are colored light-grey, bona fide GCs (c=l) are plotted in darker grey, while 
the sources that cannot be firmly classified within our scheme (c=3) are colored in the darkest grey. The box in the upper-left 
corner encloses stars and candidates whose nature requires the use of all the three parameters to be classified (V r , R, F-A flag). 
The regions of the plot corresponding to different classifications are labeled. Note that most of the clusters with R> 1.5 (Tab. 2 
and 3) do not appear in this plot for the sake of graphical clarity. 



galaxy by Racine based on ground-based, high resolu- 

tion imaging. Our morphological analysis of this object, and its 
low radial velocity, independently confirmed from Tablef^data, 
agree for a decidedly clean point-source appearance; according 
to the adopted criteria, we therefore classify it as a foreground 
star. 

4.2. Background galaxy contamination 

Since most of the known GCCs have been historically selected 
because of their fuzzy and/or extended look on photographic 
plates, it is quite natural that the major contamination source 
is represented by background galaxies, in particular those of 
spheroidal morphology. Any eye-detection survey is in fact rea- 
sonably safe with respect to low-redshift grand-design spirals, 
that can usually be confidently picked up, at least in good- 
quality images. 

Spectrophotometric information can help to segregate with 
some confidence nucleated galaxies and M 3 1 GCs in the color 
domain. In G04 we noted, for instance, that the large majority 
of spurious GCCs, eventually identified as background galax- 
ies, are redder than (V — K) > 3.0, so that by restraining tar- 
get selection to bluer objects one should in principle maximize 
the detection of genuine clusters. This guess is supported on a 
more physical basis looking at the expected photometric evolu- 
tion of early- and late-type galaxy models by Buzzoni ( 2005), 
as shown in Fig. [3] On the other hand, the contribution of 
blue galaxy contaminants (closer in color to the bulk of low- 
metallicity M 3 1 GCs) may increase as far as our GC search 
extends to fainter and more blurred targets including poorer 
candidates in the Battistini et al. image-quality classification. 



According to the updated RBC sample, the relative source 
partition of bona fide [galaxies: GCs: stars] among the entries 
with available V and K photometry and (V — K) > 3.0 is found 
to be [47: 102: 10], while the corresponding frequency blue- 
ward of the (V - K) threshold is [11: 202: 27] . Taken at their 
face value, these numbers suggest that about 2/3 of the whole 
(V - K) > 3.0 GCCs might eventually confirm to be genuine 
clusters, while our performance should raise to a nearly 85% 
for (V - K)< 3.0 candidates (see G04). 

On the other hand, one has to admit that these figures are 
at odds with the empirical evidence from the present analysis, 
as for the whole sample of 71 GCCs of Table^with firm clas- 
sification we find [galaxies: GCs: stars] = [12: 42: 17], with 
a remaining fraction of 5 Hn regions and unclassified objects 
(again, see Fig. [8] for a summary). Thus, about one in two of 
our targets eventually revealed to be a genuine globular, in spite 
of the fact that over 70% of our GCC sample is bluer than 
(V — K) > 3.0, and an a priori distribution should be expected 
such as [galaxies: GCs: stars] = [9: 57: 7] and 3 unclassified 
objects. 

These figures lead us to the following conclusions: 

(a) As far as deeper M 31 GCC surveys are carried out (and 
naturally include poorer and more "blurred" targets), the 
galaxy contamination becomes marginally more important, 
with an increasing contribution of blue spirals and star- 
forming systems. This is especially evident when we plot 
the galaxy ( V - K) distribution along the different Battistini 
et al. image-quality classification, as in Fig. [9] Fortunately, 
this bias can in principle be fully overcome as any spectro- 
scopic identification would easily discriminate these spuri- 
ous emission-line objects (note, for instance, that 9 out of 
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B-V 



Fig. 8. Apparent V - K vs. B - V color distribution for the 45 targets of Table^with available B, V, K photometry (big markers 
on the plot) compared to the corresponding RBC classification groups (small markers). Big pentagons are the bona-fi.de (c — 1 
RBC flag) M 3 1 globular clusters of this study, solid dots are background galaxies (c = 4), solid triangles mark M 3 1 Hn regions 
(c = 5), and star markers are MW stellar interlopers (c = 6) and M 31 asterisms (c = 7). The arrow indicates the reddening vector 
corresponding to E(B — V) - 0.1. The main stellar locus for stars of different spectral type (F0,F5, GO, G5 dwarfs and K0, K5, 
and MO giants from Johnson (1966) is superposed to the data, as well as the expected apparent colors vs. redshift for the Buzzoni 
(2005) template galaxy models for elliptical, spiral (type Sb) and irregular (type Im) Hubble types. Galaxy colors are tracked 
from z — to 0.5 (as labeled on the plot). The stellar locus and the galaxy models have been reddened assuming E(B - V) = 0.1 1. 
Note that the main background contaminants to M 31 GCCs (i.e. ellipticals and nucleated spirals within z < 0.2) are always 
redder than V — K = 3.0 while only star-forming spirals and irregulars (as well as MW stars) become the prevailing contaminants 
at bluer colors. 



12 recognized galaxies in our Table[0sample display some 
Balmer and/or [Oin] emission), 
(b) Between galaxy and star contaminants, the latter become 
increasingly important as far as class C-D-E Battistini et al. 
candidates are surveyed (we find 17 stellar interlopers vs. 7 
expected cases). 

4.3. Contamination by asterisms/associations 

Radial velocities and image analysis are very effective tools to 
find out the most abundant sources of contamination, i.e. back- 
ground galaxies, foreground stars and Hn regions. However, 
they can be mocked by a further and much subtler kind of con- 
taminants. Small stellar associations, very young M 31 open 
clusters (whose integrated luminosity is dominated by a few 
massive stars), and even perspective asterisms due to chance 
alignment of M 31 stars (and/or one or a few stars embed- 
ded into, or projected onto, a nebula) can easily mimics a typ- 
ical faint GCC, when observed in low spatial-resolution im- 



ages. Clearly, this class of contaminants is expected to af- 
fect faint and blue GCCs, especially those projected onto the 
disc of M 31. Cohen, Matthew & Cameron (2006) recently 
provided direct evidence of at least four such special cases 
among spectroscopically confirmed clusters, using exquisite 
spatial-resolution images obtained with adaptive optics. To ac- 
count for these "new" kind of contaminants, we introduce a 
new classification type, class c = 7, corresponding to aster- 
isms/associations ofM31 stars. 

In spite of the possible occurrence of these "fake" stellar 
aggregates, whose real occurrence needs however to be more 
firmly assessed on a statistical basis, it is clear that, in any case, 
the final word on the real nature of M 31 GCCs must come 
from the physical resolution of the composing stars, through 
very high-resolution imaging (i.e., at HST or with adaptive- 
optic ground telescopes). 

Among the clusters revealed to be asterisms by Cohen et 
al. i006), there is one, B216, we classify as a genuine GC 
in Table |2] based on its velocity and apparent morphology. 
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Fig. 9. The V-K color distribution for confirmed background 
galaxies (c = 4 flag) in the whole RBC sampl e vs. o riginal 
quality clas s according to the Battistini et al. Jl98Ci Il982l 
1 19871 1 19931) scheme. Blue (V - K < 3.0) late-type galax- 
ies (solid triangles) become an increasing contaminant source 
among class C-D cluster candidates. See text for a discussion. 

While the other asterisms in the Cohen et al. list cannot be 
independently assessed from our low-resolution imagery, both 
the BFOSC and the DSS frames for B216 appear hardly com- 
patible with the nearly empty field imaged by Cohen et al.. 
However, we carefully checked the location of this field (with 
the precious collaboration of J. Cohen) and the listed target co- 
ordinates, and we must conclude that B216 is indeed not a real 
cluster. 

4.4. Independently checked candidates 

As a duty cycle operation to maintain an updated release for the 
RBC, we periodically search the HST archive for intentional or 
serendipitous images of M 31 GCCs that can potentially reveal 
the true nature of the objects. A systematic survey of the avail- 
able material is ongoing. Here we report only a few cases in 
which a clear and indisputable confirmation can be achieved 
from a thorough inspection of the images, i.e. objects partially 
resolved into stars or obvious foreground stars or asterisms. 

- A new, clearly resolved, cluster has been identified by L.F. 
in deep F555W and F814W ACS/WFC images (see Fig. 
UOt . The cluster has no counterpart in the RBC and is 
located at a 2 ooo = 00 ;i 42 m 28.05 s , 6 20 qo = 41°33'24.5". 
According to the nomenclature adopted in G04 and G05 
we christen the newly found cluster Bologna 515 (B515). 

- the candidate B056D is clearly recognized as a genuine 
cluster on several deep (t exp up to 2370 s) ACS/WFC im- 
ages taken in different filters, by different teams. 

- NB83, classified as genuine cluster by BOO based on its ra- 
dial velocity (V r = -150 ± 14.0 km s~') is clearly recog- 
nized as a star in deep F555W and F814W images taken 
with the WFPC2. Note that the radial velocity is in the 




Fig. 10. The newly identified cluster B515 from F606W 
HST/ACS imagery. The scale and the orientation of the image 
are reported. 

range where contamination by MW stars can occur (see 
Sect. 4.1). 

- B102, classified as genuine cluster by H91 and P02 based 
on its radial velocity (V r = -235.4 + 11.7 km s~') is rec- 
ognized as an asterism formed by two stars superposed to 
a nebulosity in deep ACS/WFC images taken in different 
filters. 

- the candidate NB92 is recognized as a bright (likely fore- 
ground) star in deep ACS/WFC and in shallow WFPC2 im- 
ages taken in various filters. 

- the candidate B162, that we classified as a genuine clus- 
ter based on its radial velocity (V r = -146 + 8 km s~') 
and on its extendedness (see Tab. 3), is clearly recognized 
as a genuine cluster also on several deep (t exp up to 2370 
s) ACS/WFC images taken in different filters, by different 
teams. 

- the candidate G137, that we recognized as an HII region 
from its spectrum, appears as a bright point-source sur- 
rounded by an asymmetric nebula on deep ACS/WFC im- 
ages, thus confirming our spectroscopic classification. 

- B118 7 , NB99, NB100, NB106, classified as stars in the 
RBC are recognized as stars also in deep ACS/WFC im- 
ages, thus confirming the existing classification on a much 
firmer basis. 

In summary, we identified a new cluster (B515), one can- 
didate is recognized as genuine clusters (B056D), two ob- 
jects previously believed to be genuine clusters has been re- 
classified as foreground star and asterism (NB83 and B102, re- 
spectively), one candidate has been firmly classified as a star 
(NB92), and the existing classification of seven other objects 
(Bl 18, B162, B137, NB99, NB100 and NB106) has been fully 
confirmed. 

5. An updated sample of confirmed M31 GCs 

As a result of our spectroscopic and imaging survey of Tabled 
candidates, we have provided 42 newly confirmed bona-fi.de 

7 See Updates & Revisions http://www.bo.astro.it/M31/ 
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Fig. 11. Distribution of confirmed M31 clusters in the reddening-corrected J - K vs. B - K plane. E(B - V) = 0.1 1 is assumed 
for all M31 clusters. Small filled circles are previously confirmed clusters, pentagon markers are the newly confirmed clusters 
studied in the present paper. SSP model sequences by Buzzoni ( 1989), for fixed ages (i.e. 15 Gyr for the upper sequence and 2 
Gyr for the lower sequence) and varying metallicity (from [Fe/H] = -2 to a solar value, at equal steps of 0.5 dex, see big solid 
dots in the sense of increasing J — K color along the sequences) are superposed. A Salpeter IMF slope and a red horizontal branch 
morphology is assumed in the models. 



M 3 1 clusters, while a total of 34 GCCs should in fact be com- 
prised among background galaxies, foreground stars or Hn re- 
gions. Our study increases the total number of confirmed M 31 
GCCs from 337 to 368 and the number of confirmed GCs hav- 
ing a radial velocity estimate from 313 to 349. 

While a thorough analysis of the integrated properties (in- 
cluding Lick indices) of the newly confirmed clusters is de- 
manded to Paper II of this series (see Sect. 1), a first glance 
to their composing stellar populations is provided in Fig. 
where we compare the new GCs with the a lready confirmed 
ones (from G04) and with the Buzzoni dl989l) theoretical mod- 
els for Simple Stellar Populations (SSPs), in_the reddening- 
corrected J-K vs. B-K plane (see James et al. (2005) for a thor- 
ough discussion of this diagnostic plane). It is clear that most 
of the new clusters have the typical integrated colors of classi- 
cal old globulars, and they appear to span the whole metallicity 
range covered by previously known M 31 GCs. 

In terms of GC luminosity function, our survey provides 
a strong contribution in the 17.0 < V < 18.0 range, where 
we add 30 new clusters to the 120 previously known objects 
(see Fig. I12t . The 42 newly confirmed clusters are distributed 



along all of the Battistini at al. classes from A to D, but we 
remarkably increased (+75%) the surveyed fraction of class D 
candidates. 

It is important to consider that only 172 candidates of class 
D and E have been scrutinized to date. Since 67 of them turn 
out to be genuine clusters (-35%) and ~500 candidates of 
these classes still need to be confirmed, one might conclude 
that over 100 genuine clusters are still hidden in this harvest 
of "intermediate/low-quality" targets. Hence, large surveys, as 
the one presented in this paper, are badly needed to eventually 
reach a truly complete sample of M 31 clusters. 

5.1. A basic application: the M 31 mass estimate 

A preliminary application of the basic kinematical properties 
of our enlarged sample provides results in excellent agreement 
with the more specific analysis by P02. The mean systemic ve- 
locity of the M 3 1 GC system (after a 2-cr clipping procedure) 
is (V C c) = -296 ± 12 km s" 1 , the median is -297 ± 14 km s -1 . 
The overall velocity dispersion is <x = 158 + 10 km s _1 . 



8 Specifically, the quality-class distribution results [A: B: C: D: 
E] = [13: 10:4: 15:0]. 
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Fig. 12. The V luminosity function of confirmed M 31 clusters. 
Continuous and dashed lines show the distributions before and 
after the present survey, respectively. Most of the new addi- 
tions (15 out of 42) are for class D candidates according to the 
Battistini et al. classification scheme. 

In Fig. El we show the distribution of the M31 -centric ve- 
locities of confirmed GCs as a function of the projected dis- 
tance from the center along the major axis of the galaxy (X, 
see G04, and references therein). To obtain a more easily read- 
able plot we limited to the range -30 kpc < X < +30 kpc, 
while the outermost cluster in our catalog (B514, see GOO) lies 
at X = 59.6 kpc. 9 

Th e well known rotation pattern of M3 1 GCs (see van den 
Bergh (2000), P02, and references therein) is clearly visible in 
Fig- El with a flattening occurring for \X\ > 7 kpc. Averaging 
and 2-cr clipping velocities as a function of X in boxes 4 kpc 
wide, we obtain an amplitude of the overall rotation pattern of 
134 + 15 km s _1 , again in very good agreement with P02. 

As an example of the possible applications of the newly 
obtained large and homogeneous database of radial velocities 
of M3 1 GCs, we obtained a simple estimate of the mass of M3 1 
within R 60 Kpc, using th e Projected Mass Estimator (PME, 
Bahc all & Tremaine dl98 lh . Adopting the version by Heisler et 
al. {1985) of the PME 
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and assuming an isotropic velocity distribution (C = 32), P02 
obtained for M 31 a total mass M tot = 4.1 +0.1 x 10 11 M Q , 
by using 319 dynamical probes (GCs) out to a radius of 27 
kpc from the galaxy center. Under the same assumptions, from 
our enlarged sample of 349 GCs out to ~60 kpc from the 
galaxy center (projected distance, R p ), we obtain M tot = 4.4 + 
0.2 x 10 11 M . It has to be noted that the uncertainties on the 
actual isotropy degree of the underlying velocity distribution 
contributes an additional factor ~2 to the quoted uncertainties. 
Note also that the estimate is unchanged if we exclude from the 
samp le the presumably young clusters identified in Fusi Pecci 
et al. ( 2005) as belonging to the thin disc of M31. The obtained 
value is well within the range spanned by the most recent es- 
timates of the mass of M 31, as listed by Evans & Wilkinson 

9 For M 31 we adopt a distance modulus (m - M) = 27.47, from 
McConnachie et al. J2005), and E(B-V)= 0. 1 1, as in G04. This corre- 
sponds to a distance D = 783 kpc. 
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Fig. 13. Radial velocity of the M3 1 globular clusters (corrected 
for the systemic velocity of M31) vs. the projected distance 
along the major axis (X). Pentagon markers represent the newly 
confirmed clusters studied in the present paper. The continuous 
line is a fit to the 2-cr clipped mean velocity as a function of 
major-axis projected distance, computed on 4 kpc wide boxes 
shifted by 1 kpc steps in X. 2-cr contours are plotted as dotted 
lines. 



(2000) in their Table 6. The agreement with previous estimates 
obtained using GCs as tracers (Federici et al dl993l) . P02) is 
also quite good. 

The Bahcall & Tremaine 



1981) PME was originally con- 
ceived for test particles to probe a central point-mass gravita- 
tional source. As an illustrative exercise, we can suitably ap- 
proach this condition by restraining our analysis to the 14 most 
distant clusters in our sample, with R p > 20.0 kpc, and main- 
tain the isotropy hypothesis (i.e. C = 16, in eq^above)- With 
these constraints we obtain M tot = 4.3 - 7.0 x 10 11 M Q , where 
the reported range has bee n obtained by a jackknife resampling 
technique (Lupton ( 1993)). This estimate is in good agreement 
with the results by Federici et al. (1993) and with the recent 
independent estimates by Carignan et al. (2006) and Chapman 
et al. ( I2006I) . 

The above estimates rely on methods that assume that the 
adopted tracers follow the mass density profile of the probed 
potential, that is not the case for GC s syste ms, in general. To 
overcome this problem Evans et al. (2003) introduced a new 
mass estimator that doesn't require any coupling between the 
distribution of the tracers and the underlying mass distribution. 
As an example of application, Evans et al. (hereafter E03) pro- 
vide an estimate of the mass of M3 1 based on GCs, obtaining 
M — 1.2 x 10 12 M . With the same assumptions as Evans et al. 
we obtain M rot = 2.9 x 10 11 M for the rotational component 
(to compare with M rot = 3.0 x 10 11 M found by E03) and 
M pres = 2.1 x 10 12 M for the pressure component (to compare 
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with M pres = 0.9 x 10 12 M found by E03). The total mass is 
M = 2.4 x 10 12 Mq, a factor 2 larger than the E03 estimate. Part 
of this difference is due to the larger value of the M3 1 distance 
adopted here with respect to E03. Our experiments, however, 
indicates that the results from this method are quite sensitive to 
the the way in which the rotational and pressure components 
are disentangled: the very simple rotation pattern adopted here 
(see Fig. ll3t introduces a large uncertainty in our result. Once 
the above factors are taken into account the agreement with E03 
is satisfying, at least in this preliminary stage of the analysis. 

A detailed analysis of the kinematics of the M 3 1 GC sys- 
tem and the galaxy mass profile is deferred to the completion 
of ou r remo te-cluster search, currently in progress (see Galleti 
et al. d2005l) l. 

6. Summary and conclusions 

We have presented the first results of a large spectroscopic and 
imaging survey of candidate clusters in M 31. The survey al- 
lowed us to confidently classify 76 candidates whose nature 
was previously unknown: 42 of these resulted new M 3 1 GCs, 
while 12 have been recognized as background galaxies, 2 are 
M 31 Hn regions, while the remaining 17 objects are fore- 
ground stars and 3 unclassified objects (possibly M 3 1 clusters 
or foreground stars). An estimate of the radial velocity has been 
obtained for all the 42 newly recognized clusters as well as for 
an additional sample of 55 M 3 1 GCs previously confirmed by 
other authors and two controversial objects. The various set of 
radial velocities for M 31 GCs available in the literature has 
been reported to the same scale, multiple measures have been 
averaged (see Sect. 3.2 for details), and a final merged catalog 
has been produced (see Tab. 1, Online material). The present 
analysis has increased the sample of confirmed M 31 cluster 
from 337 to 369 members, and the number of confirmed GCs 
with a radial velocity estimate is increased from 313 to 349. 

While the main basis for the classification work was pro- 
vided by radial velocities, we have also implemented a method 
to distinguish between point and extended sources on low- 
resolution imaging that allows - in many cases - to disentangle 
genuine M 3 1 clusters and MW stars when radial velocity alone 
leads to controversial conclusions. We also provide a safe clas- 
sification for few candidates not included in our survey based 
on the inspection of publicly released high spatial resolution 
images from the HST general archive. 

6.1. The Revised Bologna Catalog V2.0 

All the present observational material has been consistently 
implemented to update the RBC, now available on line in 
its latest V2.0 release. Future minor updates of the cat- 
alog will be described and commented in the RBC web 
page (http://www.bo.astro.it/M31/) where all the cata- 
log database can be retrieved as ASCII files. 
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